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Abstract 

Over the last few years numerous papers concerning analog mod- 
els for gravity have been published. It was shown that the dynamical 
equation of several systems (e.g. Bose-Einstein condensates with a 
sink or a vortex) have the same wave equation as light in a curved- 
space (e.g. black holes). In the last few months several papers were 
released which deal with simulations of the universe. 

In this article the de-Sitter universe will be compared with a freely 
expanding three-dimensional spherical Bose-Einstein condensate. Ini- 
tially the condensate is in a harmonic trap, which suddenly will be 
switched off. At the same time a small perturbation will be injected 
in the center of the condensate cloud. 

The motion of the perturbation in the expanding condensate will 
be discussed, and after some transformations the similarity to an ex- 
panding universe will be shown. 
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1 Ultra-cold dilute gases 



The Hamiltonian for a ultra-cold gas — consisting of Nq molecules — in a 
time-dependent external trap is [I] 



H = / df^(t,r) 



-^v 2 + v ext {tA if(t,f) 

+ 2 ^dr' #(t, r) $f(t, f)' V(f- r') &(t, ?) &(t, r) . 

The terms in the first line express the kinetic energy and the external po- 
tential, while the second line describes the interaction part, which only takes 
into account two-particle interactions. It is well known that the evolution of 
a field operator is given by the Heisenberg-relation: 

ih^(t,r) = $(t,r),H]. (2) 

Restricting the discussion to the case of zero temperature, the ground state 
of the condensate can be characterized by the Hartree-Fock ansatz 

|$) = |$) ®...<gi |$), (3) 

where the field operator iff is replaced by Nq identical macroscopic wave- 
functions <E>, and in that regime the interaction potential can be replaced by 
a so-called pseudo potential. The model is a hard-core potential, which is 
completely described by the scattering length a. Here a depends on serval 
factors — for example the external field, what kind of atoms are be used and 
so on — but it can experimentally measured. Together with the normali- 
sation condition, the ground state of the macroscopic wave-function has to 
fulfill the Gross-Pitaevskii equation: 

iTi = ("|^V 2 + V ext {t,r) + N og mr}\ 2 - ^ $(t,f) , (4) 

where g = a is the coupling constant. 



2 Sound waves in the condensate 

Small laser- injected density modulations in the condensate cloud had first 
been measured by Andrews et al. 0. If the perturbation is small enough 
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the macroscopic wave-function remains in the ground state. The dynamics of 
the perturbation, in a time-dependent external trap, is given by the Gross- 
Pitaevskii equation. In the eikonal approximation the density p(t, r), and the 
phase 9(t, r), can be separated from each other in the following way [HHHE]: 



$ = y/p(t, r) exp{z 9{t, r)}. 



(5) 



A small perturbation in the cloud can then be characterized by a small vari- 
ation in the density p = p + ep\, and the phase 9 = 9 + e9\. Using the 
Gross-Pitaevskii equation leads to 
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-V(p V9 1 + Pl V9 



Here the so called Thomas-Fermi equation has been used. In the ground 
state, the kinetic energy of the condensate cloud is small compared to the 
potential and interaction energies. 

These two equations can easily be rewritten as a single equation for the 
phase of the perturbation El E] : 



-0i - d t {vV9 x ) - V{v9\) + V ((c 2 - v 2 )V9 1 ) = 



(7) 



Introducing an effective metric g^, the differential equation for the pertur- 
bation of the phase has the same structure as the wave-equation [SI El El El : 



for a massless particle in a curved-space-time. The effective metric is 



(8) 
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Here c, the speed of sound and v, the background velocity are given by 



(9) 



gpo 



in 



(10) 
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and 

v = -ve (ii) 

m 

respectively. 

Considering the last four equations, we see that, once the speed of sound 
and the background velocity are estimated, the dynamics for that config- 
urations is fully given by wave equation (JSJ). We now use these results to 
investigate a time-dependent expanding condensate. 



3 Expanding condensate cloud 

Initially the condensate is in the ground state. This is achieved by a static 
three-dimensional harmonic trap. Suddenly, at time t = the external trap 
is switched off. Immediately, the condensate cloud will expand. 
In [7j it is shown, that all small portions of the gas expand along trajectories 

r{t) = z{r,t = 0)-b(t). (12) 

Consider an infinitesimal shell at the radius r = r\ and at the time t = 0. 
After a time interval t the shell radius will be increased to r(t) = z(ri,t = 
0) • b(t). It is convenient to choose r(t, 0) = z(r,t = 0) corresponding the 
initail configuration, therefore b(0) has to be unity. In Fig.(£Q) the expansion 
of a single shell is illustrated. 
Then, the radial background velocity is 

v T = ^r. (13) 

For the speed of sound the proportionality to the density can be used, because 
the expansion for the volume, with respect to time, is V(t) = b(t) 3 ■ V(0). 
Then the speed of sound is given by 

*r) = ^- (14) 

Next, the results from this section and those from the previous section 
will be put together. 
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condensate at time t = 

condensate at time t = t\ 

r(ti) = 6(tiMr,Q)____ 




Figure 1: For times t < the condensate is in the ground state. At time t = 
the external trap will suddenly be switched off and the expansion starts. 

4 de— Sitter universe 

In spherical coordinates the line element for the metric of Eq.© is jH] 

ds 2 = c(r,t) [(-c(r,t) 2 - v(r,t) 2 ) dt 2 + 2v(r,t) dr dt + dr 2 + r 2 dtt 2 ] , 

(15) 

where Q is the angular-space element. 

Rewriting the line element in the coordinates chosen for the trajectories leads 
to 

ds 2 = c (r) [-co(t) 2 b{ty 9/2 dt 2 + b(t) l/2 dz 2 + b{tfl 2 dz 2 dQ 2 ] . (16) 
Here co(r) = c(r, t = 0) is the speed of sound at the time t = 0. 

To compare this line element with the de-Sitter line element one last 
transformation is necessary, which eliminates the time-dependence of dt 2 . 
The appropriate transformation for that is given by 

dr = b(ty 9/4 dt , (17) 

and the resulting line element is 

ds 2 = c(r) [-c(r) dr 2 + b(t) 1/2 r 2 dVl 2 ] . (18) 
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For b(t) = (9 H t + 1) 4 / 9 — where H is the well-known Hubble constant — the 
perturbation in the condensate obeys the same dynamic as the predictions 
for light in the de-Sitter universe. 

To obtain these results, two restrictions have to be made: 

• The first is a restriction to time, because the Thomas-Fermi approxi- 
mation — used in order to neglect the kinetic energy — becomes invalid 
for late times: During the expansion the potential and interaction en- 
ergies will be transformed into kinetic energy. 

• Second, the derived line element is, because of its dependence on the 
speed of sound, still space-dependent. But, for perturbations around 
the point r = 0, the speed of sound is approximately constant. In this 
approximation the line element also space- independent. 

Other approaches to the use of Bose-Einstein condensates to simulate the 
Friedmann-Robertson- Walker universe can be found in references |9~| HUj. 

5 Discussion 

The fundamental questions are: "Is this system realisable in the labaratory?" 
and "What does it teach us?" 

In contrast to other models considered in the literature, this is one where 
the experimental techniques are already known. Such experiments, called 
time of flight measurements of expanding condensates have been done in the 
laboratory. One main issue causing significant troubles is the instability. 
This problem could be solved using another configuration, e.g. using a time- 
dependent scattering length a(t), instead of an expanding condensate cloud. 
At least, we would have a tool to simulate different scenarios for the expansion 
of the universe. 
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